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Critical field in a model with local pairs 
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We analyze the role of Zeeman and orbital pair breaking mechanisms in models appropriate 
for short coherence length superconductors. In particular, we investigate the attractive Hubbard 
and the pair hoping models. The orbital pair breaking mechanism dominates in the majority of 
models with s-wave and d-wave superconducting order parameters. On the other, the repulsive 
pair hopping interaction leads to 77-type pairing, that is stable against the orbital pair breaking. 
External magnetic field reduces this type of pairing predominantly due to the Zeeman coupling. 
According to the recent experiments this mechanism is responsible for closing of the pseudogap. 
Moreover, the temperature dependence of the gap closing field in r^-phase fits the experimental 
data very well. We discuss whether the preformed pairs in the 77-phase could be responsible for the 
pseudogap phenomenon. 



I. INTRODUCTION 

Almost twenty years after discovery of high- 
temperature superconductors (HTSCs) the mechanism 
responsible for their unusual properties remains unclear. 
The complex phase diagram of HTSCs suggests that 
there may be no single mechanism that dominates over 
the entire doping range. In particular, the normal-state 
properties in underdoped and overdoped regimes are dif- 
ferent. Highly overdoped compounds in the normal state 
exhibit Fermi liquid behavior, whereas the superconduct- 
ing state may be described within a weak-coupling BCS 
theory T]. On the other hand, in the underdoped regime 
the HTSCs exhibit unconventional features. The most 
remarkable of them are the extremely short coherence 
length and a pseudogap that opens in the normal state. 
The presence of the pseudogap has been confirmed with 
the help of various experimental techniques like: angle- 
resolved photoemission Jp, 0, , intrinsic tunneling spec- 
troscopy [101, NMR UHl, infrared and transport 
measurements. Although, there is no complete the- 
oretical description of the pseudogap, one usually con- 
siders this phase as a precursor of the superconductivity. 
According to this hypothesis formation of Cooper pairs 
starts at temperature T*, higher than the superconduct- 
ing transition temperature Tc- Then, at Tc, these pre- 
formed pairs undergo Bose-Einstcin condensation. 

This hypothesis seems to be supported by recent ob- 
servations of the vortex-like Nernst signal above Tc ^3 
that evolves smoothly into the analogous signal below the 
superconducting phase transition 12]. The Meissner ef- 
fect does not occur in the pseudogap phase due to strong 
phase fluctuations rather than the vanishing of the su- 
perfluid density. Therefore, theoretical description of the 
suppression of the Meissner effect requires an approach 
beyond the mean-field level. Despite the absence of the 
Meissner effect above Tc, one can observe inhomogeneous 



magnetic domains that are interpreted as precursors to 
the Meissner state [T^ . 

The short coherence length indicates that the pair- 
ing takes place in the real space, leading to boson-like 
objects. A few models are commonly used to describe 
systems with the loc al p airs. Namely, the attractive 
Hubbard (AH) model |l4j . fermion-bosonlTBI and purely 
bosonic models;16i]. as well as the Penson-Kolb (PK) 
model0|, i.e., the tight-binging model with local pair 
hopping. These models should be considered as effective 
approaches which do not explain the microscopic origin 
of the pairing interaction. 

Another unusual property of HTSCs is related to their 
behavior in the external magnetic field. In particular, 
temperature dependence of the upper critical field Hc2 
has a positive curvaturejlSl [l9| in contradistinction to 
classical superconductors, where a negative curvature is 
observed. Moreover, Hc2 does not saturate even at gen- 
uinely low temperature. Recent experiments |2Cll| show 
that also the pseudogap is destroyed by sufficiently high 
magnetic field, Hpg. Although, the temperature depen- 
dence of Hpg has a negative curvature, it significantly 
differs from the predictions of the standard Helfand- 
Werthammer theory [22. Namely, Hpg{T) has a large 
slope at temperatures close to T* and saturates already 
at T ~ 0.7r*. These features may assist in verification 
of the preformed Cooper pairs hypothesis and, more gen- 
erally, in choosing the most appropriate model of HTSC. 



II. MODEL 

In the present paper we show that opening of the pseu- 
dogap and its dependence on the magnetic field can be 
described within a model with local pair hopping. Our 
starting point is the two-dimensional (2D) Penson-Kolb 
model with the Hamiltonian given by: 
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_lj^e2'*--4,c,V,xc,T. (1) 

Here, c^^J creates (annihilates) an electron with spin a at 
site i, tij is the single electron hopping integral between 
sites i and j, //is the chemical potential and J is the 
nearest neighbor pair hopping interactions. The external 
magnetic field perpendicular to the lattice Hz shifts the 
energy levels by guBH^a (g is the gyromagnetic ratio 
and fj,B is the Bohr magneton) and modifies the hopping 
terms. The single electron hopping integral acquires the 
Peierls factor 

,.,^±Qa.,1 (2, 

whereas the phase factor in the pair hopping term is twice 
larger. 

The Penson-Kolb model can be derived from a gen- 
eral microscopic tight -binding Hamiltonian j22j| . where 
the Coulomb repulsion may lead to the pair hopping in- 
teraction. In such a case J is negative (repulsive Penson- 
Kolb model). However, treating the Penson-Kolb model 
as of a phenomenological nature, we assume J to be an 
effective parameter, that can be negative as well as posi- 
tive. It can be understood as a result of renormalization 
originating, e.g., from electron-phonon coupling 23]. For 
a nonzero single electron hopping integral J — > — J is not 
a symmetry of the PK model j2J|. However, supercon- 
ducting correlations occur in the Penson-Kolb model for 
attractive pair hopping interaction (J > 0) as well as for 
the repulsive one ( J < 0), provided that the pair hopping 
is large enough. The latter case is usually referred to as 
ry-type pairing. Then, the total momentum of the paired 
electrons is Q = (tt, tt) and the phase of superconducting 
order parameter alters from one site to the neighboring 
one. It has been shown that there is a flux quantization 
and Meissner effect in this state. 25] Superconductivity 
survives also in the presence of on-site Coulomb repul- 
sion (Penson-Kolb-Hubbard model), provided that this 
interaction is not too large p^. 

A. Density of states 

At the mean-field level, for J > one obtains an 
isotropic superconducting gap, identical to that obtained 
for AH model. On the other hand, in the case of 77-type 
pairing (J < 0), the density of states is finite for arbi- 
trary energy. However, the density of states at the Fermi 
level may significantly be suppressed for some dopings. 
In the simplest case of the nearest neighbor hopping the 
density of states in the 77-phase is of the form: 

P(w) = ^ (1 - 1) Po(^-m) + ^ (1 + 1) Po(^+A)- (3) 

Here, fi = fi'^ + AlJA]^ , A = (-1)* {c^lC,^) is the rj- 
phase order parameter and po is the density of states for 
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FIG. 1: Density of states for the Penson-Kolb model in the 
77-phase. We have used t' = —0.25t, J = —t and fi = — 0.75t. 
The curves from the topmost to the lowest correspond to the 
values of the order parameter A = 0, 0.1, 1. The dashed 
line indicates the Fermi level. 

^ — A — 0. One can see from Eq. Q that the quasi- 
particle poles split when A becomes finite. Therefore, 
a local minimum in the density of states may occur at 
the Fermi surface. Despite the presence of this minimum 
the density of states at the Fermi level remains finite 
provided that p, is small when compared to the band 
width. Inclusion of the next nearest neighbor hopping t' 
leads to a more complicated expression for the density 
of states. However, the structure of p in the ry-phase re- 
mains unchanged. Fig. 1. shows the density of states 
calculated for t' ^ and different values of the Ty-phase 
order parameter. Gradual decreasing of p at the Fermi 
level resembles opening of the pseudogap in HTSCs. 

Another feature that could speak in favor of this in- 
terpretation is anisotropy of the gap Q. More precisely, 
for f ^ the magnitude of splitting of the quasiparticle 
peaks depends on the direction in the Brillouin zone. The 
splitting of the spectral functions is presented in Fig. 2. 
As we consider isotropic order parameter the splitting is 
finite everywhere at the Fermi level, in contradistinction 
to a purely rf-wave gap. However, this drawback may be 
removed when considering a nonlocal pairing. 

B. Response to magnetic field 

In contradistinction to the AH model, the exter- 
nal magnetic field explicitly enters the term responsi- 
ble for superconductivity, i.e., the pair hopping inter- 
action. Therefore, the differences between AH and PK 
models may show up in the electromagnetic properties 
|27j ]. Here, we investigate the temperature dependence 
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FIG. 2: The spectral functions in various points at the Fermi 
surface. We have used the same model parameters as in Fig. 
1. The position of points A and B is depicted in the inset. 
The continuous curve corresponds to the case A = 0, whereas 
the dashed lines have been obtained for A = 0.5. 

of critical field i?crit- It is defined as the highest magni- 
tude of the magnetic field, for that there exists a non-zero 
solution for the superconducting order parameter: 

A. = (c^iQ|). (4) 

As we carry out calculations at the mean-field level it 
is impossible to determine the phase coherence of the 
Cooper pairs. Therefore, the physical interpretation of 
the critical field is not unique. In the overdoped regime 
of cuprates it can directly correspond to the upper crit- 
ical field. On the other hand, for underdoped systems, 
and within the precursor scenario of pseudogap, it can 
be interpreted as ^/pg, i.e., the field at which the inco- 
herent pairs appear. One can argue that the PK model 
can also be used as an effective model of other super- 
conductors that are characterized by a short coherence 
length: barium bismuthates, fuUerides, Chevrel phases, 
organic superconductors and heavy fermion systems, 
In these cases i?ciit should be considered as ifc2- 
For the sake of simplicity we define: 

i 

where the prime means that the summation is carried out 
over the nearest neighbors of site i. Then, the mean-field 
Hamiltonian takes on the following form: 

i,j.(7 i(7 



At the mean-field level the only difference between PK 
and AH models is the presence of in Eq. © instead 
of Ai. Therefore, in order to calculate the critical field 
one can follow an approach, that has previously been 
developed for the lattice gas with on-site attraction . 
Then, one ends up with the lattice version of the Gor'kov 
equations: 

A.. = ^E AjG(z,j,w„)G(i,j,-w„). (7) 

Here, G{i,j,LUn) is the one-electron Green's function in 
the presence of a uniform and static magnetic field. w„ 
denotes the fermionic Matsubara frequency. With the 
help of Eq.© one can eliminate A^ from the Gor'kov 
equations. Then, -ffcrit(T') can be calculated from 

= ^ E e2**'^A,G(z,;,c^„)G(z,;,-a;„) (8) 

or 

= ^ E (-l)^+'e2^*'^A,G(z,/,c.„)G(z,;,-a;„). 

(9) 

Eqs. lO and © are equivalent since (— 1)-'+' ~ — 1, 
for the neighboring sites j and I. However, it is more 
convenient to use the first/second of them for attrac- 
tive/repulsive pair hopping interaction. 

In the following we consider only the nearest-neighbor 
one-particle hopping integral t and use the Landau gauge 
A = Hz (0,a;,0). Then, the Harper equation 

9 iPx,Py, a; + 1) + 2 cos {hx ~ pyo) g {px,Py,x) 
+g{Px,Pv,x- 1) = r'^E{px,py)g{px,py,x), (10) 

determines eigenvalues E {px,Py) of the one-particle hop- 
ping term. The corresponding eigenstates are enumer- 
ated by PxiPy and are of the form: 

Ux,y {Px,Py) = e^Pyy°-g{Px,Py,x) . (11) 

Here, a:, y are integers which enumerate the lattice sites 
in X and y directions, whereas h/(2'K) is a ratio of the 
flux through a lattice cell to one flux quantum. We refer 
to Ref. ^] for the details. 

The one-electron Green's function can be expressed 
with the help of eigenvalues and eigenstates of the 
normal-state Hamiltonian. Then, the summation over 
Matsubara frequencies in Eqs. lO and ^ can explic- 
itly be carried out. In the Landau gauge the presence of 
magnetic field does not change the plane-wave behavior 
in y-direction [see Eq. ifTTjl ]. Therefore, the supercon- 
ducting order parameter depends only on x and pairing 
of electrons takes place for the same y-components of 
their momenta, as in the absence of magnetic field, i.e., 
{py, —py) for J > and {py, tt — py) for J < 0. Taking 



these features into account one can rewrite the Gor'kov 
equations for the attractive: 

X [2cos(2/ix) g{px,Py,x)g {k^,-py,x) 

+ 9 {Px,Pv,x + l)g{k^, ~py,x + 1) 

+ 9{Px,Py,x-l)g{h,-Pv,x-l)] 

'x g(.Px,Pv,x')g{h,-Py,x') , (12) 

as well as for the repulsive pair hopping interaction: 

n^:' = ^^Xl^^r X{Px,Py,h,T^ -Py) 

Py iPx J 

X [2cos(2/ia;) g {px,Py, x) g (kx,Tr - Py,x) 

- g iPx,Py,X + l)g (fcx, Tr-py,x + l) 

- g{Px,Py,X-l)g{h,Tr-Py,X-l)] 

X g{Px,Py,x')g{h,Tr-Py,x') . (13) 

Here Hx = Ax^y{—iy and the Cooper pair susceptibility 
xiPx^Py, kx,ky) has the following form: 



X{Px,Py] kx, ky) 



tanh 



tanh 



E {px,Py) - ^J.- g^isH^ 
2kBT 

E {kx, ky) - 11 + gusH^ 



2kBT 

X [2{E{px,Py)+E{kx,ky) -2n)y . 

(14) 

The above equations determine the strength of mag- 
netic field at which the local pairing disappears. We have 
carried out calculations for 150 x 150 cluster with peri- 
odic boundary conditions (be) along the y axis. As the 
Landau gauge breaks the translation invariance along x 
axis we have used fixed be in this direction. Our pre- 
vious calculations indicate that such a size of cluster is 
sufhcient to obtain convergent results. [28j 

Fig. 3 shows the temperature dependence of i/crit ob- 
tained for attractive and for the repulsive pair hopping 
interaction. These results are compared with Hcrit{T) 
calculated from the 2D AH model|2^ with U = —t. We 
have adjusted the strength of the pair hopping interac- 
tion to obtain the same critical temperature in the ab- 
sence of magnetic field. For J > ffcrit(r) in PK model 
is very close to that of AH model. It means that in 
the case of s-wave pairing the Peierls factor in the pair 
hopping term leads only to a small decrease of supercon- 
ducting correlations. However, the temperature depen- 
dence of Hcrit in the ry-state differs qualitatively from 
the s-wave case. Namely, -ffcrit {T) has a very large slope 
for a weak magnetic field and saturates already at rela- 
tively high temperature. Such a behavior of the critical 
field resembles Hpg{T), that has recently been observed 
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FIG. 3: Temperature dependence of ffcrit for t' = and /x = 0. 
Continuous curve has been obtained for the PK model with 
J — —1.56t (ry-type pairing). The dashed line corresponds 
the attractive pair hopping interaction J = 0.5t. The dotted 
line shows the critical field in the AH model with U = t. 
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FIG. 4: Fit of the theoretical results to the experimental 
data for the pseudogap closing field -H^P9(T).[23| The contin- 
uous line represents i/crit(T) calculated for the repulsive pair 
hopping interaction. We have used the same model parame- 
ters as for the r;-pairing in Fig. 3. 



in Bi2Sr2CaCu208+y.^i In Fig. 4 we compare our re- 
sults and the experimental data. 

Within the Helfand-Werthammer theory, the temper- 
ature dependence of critical field is predominantly de- 
termined by the diamagnetic pair breaking mechanism. 
The Zeeman coupling becomes important only for suf- 
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FIG. 5: Hciit{T) calculated for attractive and repulsive PK 
model with and without the Zeeman coupling. We have used 
the same model parameters as in Fig. 3. 

ficiently strong magnetic field. This feature fiolds also 
in the case of the lattice gas.|2^ In order to investigate 
the role of the Zeeman and orbital contributions in PK 
model, we have repeated our calculations in the absence 
of the Zeeman term. The resulting iJcrit {T) is shown in 
Fig. 5. 

In contradistinction to the s-wave superconductivity, 
the diamagnetic pair breaking is of minor importance 
in the case of 77-pairing. This feature is responsible for 
extremely high values of Hcrit in the absence of Zee- 
man term. Experimental investigations^Oj show that the 
pseudogap closing field scales linearly with T*. In Ref. 
|20l | the value of the scaling factor has been interpreted in 
favor of the Zeeman coupling as a mechanism that closes 
the pseudogap. This pair breaking mechanism dominates 
also in the case of 77-pairing. 

III. SUMMARY 

To summarize, we have investigated the PK model 
with attractive as well as repulsive pair hopping interac- 
tion. We have shown that the repulsive pair hopping term 



may lead to the occurrence of local minimum in the den- 
sity of states, that is characteristic for pseudogap phase 
of underdoped cuprates. It originates from the spitting of 
the quasiparticle peaks. Despite the on-site pairing the 
magnitude of the splitting is a direction-dependent quan- 
tity, provided that t' ^ 0. Anisotropy of the pseudogap is 
observed in ARPES experiments. Q We have also calcu- 
lated the temperature dependence of iJcrit , defined as the 
highest magnetic field for which there exists a non-zero 
solution for the order parameter. We have found that 
in the case of 77-type pairing i?crit (T) reproduces the ex- 
perimental data for the pseudogap closing field. These 
features do not occur for attractive pair hopping interac- 
tion. In this case the gap structure as well as iJcrit(T) 
are similar to those obtained for AH model. 

Our approach to the critical field accounts both for 
Zeeman and diamagnetic pair breaking mechanisms. In 
the case of s-wave pairing inclusion of the Zeeman cou- 
pling does not lead to any essential changes in Hcnt{T). 
On the other hand, Zeeman term is of crucial importance 
for ry-pairing, whereas the diamagnetic pair breaking is 
ineffective. According to the experimental data the pseu- 
dogap is closed by the Zeeman splitting. 

As we have previously|2^ shown, other models appro- 
priate for short coherence superconductors have ground 
states (with s-wave or d-wave symmetry) that are al- 
most insensitive to the Zeeman interaction. Therefore, 
the PK model with J < is unique in that the gap is 
closed predominantly due to the Zeeman interaction. 

Collecting the features: the presence of the pseudo- 
gap, its anisotropy, Zeeman origin of -ffcrit (in agreement 
with the experimental data), the presence of flux quanti- 
zation and the Meissner effect, (consistent with the pre- 
formed Cooper pairs scenario), may lead to a tempting 
hypothesis that the pair hopping can be responsible for 
the pseudogap. However, in order to avoid the problem 
of interpretation of the critical field, it should be verified 
beyond the mean-field level, discussed in this paper. 
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